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ABSTRACT 


A  dynamic  instability  is  described  for  rolling,  finned  missiles  with 
canted  fins.  The  instability  occurs  as  an  undaaplng  of  the  angle  of  attack 
and  is  caused  by  differential  lift  from  unequal  effectiveness  of  the  windward 
and  leeward  fins.  The  instability  is  similar  to  a  Magnus  instability  in  that 
a  yawing  moment  occurs  that  causes  damping  of  negative  precession  and  undamplng 
of  positive  precession  motions.  The  equivalence  of  the  Magnus  coefficient  and 
the  differential  fin  lift  forces  is  derived,  and  the  instability  is  demonstrated 
with  computer  simulations  of  the  equations  of  motion  for  the  angle-of-attack 
convergence  of  a  reentry  vehicle.  The  computer  results  are  compared  with  a 
closed-form  solution  for  the  angle-of-attack  convergence  envelope. 
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HCH£mCLATuk£ 

constant 

slope  of  fin  lift  effectiveness  vs  angle 
of  attack 

aerodynamic  fin-induced  roll  moment 
coefficient 

aerodynamic  fin  lift  derivative 
aerodynamic  pitch  damping  derivative 

aerodynamic  Magnus  moment  derivative 

aerodynamic  normal  force  derivative 

aerodynamic  reference  diameter 
pitch  or  yaw  moment  of  inertia 
roll  moment  of  inertia 

fin  lift  force 
vehicle  mass 
aerodynamic  yaw  moment 

aerodynamic  pitch  moment 

aerodynamic  roll  moment 

roll  rate 
dynamic  pressure 
aerodynamic  reference  area 
vehicle  velocity 


PRECOOK  PACE  HANK 


NOMENCLATURE  (Continued) 


x  -  distance  of  fins  aft  of  center  of  mass 

x  -  static  margin  (distance  of  center  of 

&  pressure  aft  of  center  of  mass) 

(x,y,z)  ■  body-fixed  coordinates 

(X,Y,Z)  *  inertial  coordinates 

6  »  effective  fin  cant  angle 

6  *  fixed  fin  cant  angle  relative  to  body 

°  axis  of  symmetry 

f,  -  yaw  axis 

il  «  pitch  axis 

0  -  total  angle  of  attack  (Euler  angle) 

<  -  roll  moment  coefficient,  C.  qSd/I 

*o  x 

li  ■  inertia  ratio,  WI 

v  “  pitch  or  yaw  damping  coefficient, 

(qSd2/2iu)(-C  +  2C  1/md2) 

m  w 

q  a 

£  «*  roll  axis 

a  *  2w/up 

T  «  C5  qSx/I 

<p  ■  roll  orientation  relative  to  wind  (Euler 

angle) 

ip  -  precession  angle  (Euler  angle) 

u  -  natural  pitch  frequency 


I.  INTRODUCTION 


It  has  been  observed  in  flight  and  demonstrated  mathematically  that  a 
statically  stable  missile  at  zero  roll  rate  can  become  dynamically  unstable 
at  sufficiently  high  roll  rates  as  a  result  of  Magnus  forces.*  *  In  effect,  the 
Magnus  forces  induce  positive  damping  that,  for  sufficiently  high  roll  rates, 
overcomes  the  negative  yaw  damping  that  is  usually  present  to  some  degree. 

The  Magnus  force  is  characteristically  dependent  on  both  angle  of  attack  and 
roll  rate,  analogous  to  the  classical  Magnus  lift  on  a  spinning  cylinder  in  a 
cross  flow.  As  such,  Magnus  instabilities  have  been  observed  on  both  finned 
and  unfinned  bodies  of  revolution. 

The  addition  of  fins  to  an  axisymmetric  missile  can  alter  its  stability 

characteristics  by  Increasing  the  classical  Magnus  forces  or  by  causing  Magnus- 

4  5 

type  effects  that  have  a  different  origin  than  the  classical  Magnus  forces.  ’ 

4 

It  has  been  demonstrated  from  wind  tunnel  tests  that  body  interference  on 
a  finned  missile  at  angle  of  attack  can  cause  an  unbalance  in  fin  lift  forces 
that  produces  a  net  yawing  moment  analogous  to  the  body-induced  Magnus  moment. 
The  fin-induced  yaw  moment  can  act  in  the  same  or  opposite  direction  to  the 
body-induced  Magnus  moment,  depending  on  the  angle  of  attack  and  on  whether  or 
not  the  missile  has  reached  Its  terminal  roll  rate  dictated  by  the  fin  cant 
angle  and  freestream  velocity. 

This  paper  describes  a  mathematical  model  for  predicting  a  dynamic 
instability  similar  to  that  identified  in  Ref.  4.  This  model  assumes  that  an 
unbalance  in  fin  lift  forces  for  a  missile  with  canted  fins  arises  from 
unequal  effectiveness  of  the  windward  and  leeward  fins  when  the  missile  is 
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at  angle  of  attack.  The  equivalent  y-  damping  and  Magnus  coefficient  are 
obtained  as  functions  of  the  fin-induced  roll  moment  and  the  relative 
effectiveness  of  the  windward  and  leeward  fins.  A  closed-form  solution  is 
obtained  for  the  angle-of-attack  convergence  envelope  of  a  rolling  reentry 
vehicle  with  canted  spin  fins.  This  is  an  extension  of  previous  work,  in 
which  a  solution  was  obtained  for  the  angle-of-attack  convergence  envelope 
of  a  nonfinnned  reentry  vehicle  with  roll  acceleration  and  pltch/yaw  damping. 


II.  EULER  ANGLE  COORDINATES 


The  vehicle  rotational  motion  is  described  in  terms  of  the  Euler 

angles  ip,  0,  6,  which  describe  the  position  of  a  set  of  body-fixed  axes 

x,  y,  s  relative  to  an  inertial  frame  X,  2.  that  translates  with  the 

vehicle,  as  shown  in  Fig.  1.  The  axes  £,  n>  C»  are  axes  of  roll,  pitch, 

and  yaw,  respectively,  relative  to  the  plane  of  total  angle  of  attack  They 

precess  about  the  velocity  vector  with  angular  rate  ip.  The  Euler  angle  6 

is  the  total  angle  of  attack,  and  the  angle  <j>  is  the  roll  angle  relative 

to  the  wind.  The  roll  rate  p  is  then  the  roll  rate  relative  to  the  wind  $ 

« 

plus  the  component  of  precession  ip  cos  6  along  the  roll  axis;  i.e., 

p  -  $  +  Ip  cos  0  (1) 

If  the  principal  moments  of  inertia  about  the  £,  n,  C  axes  are  1^,  I,  I, 

respectively,  and  the  aerodynamic  moments  about  these  axes  are  M, ,  M  ,  M  , 

s  n  s 

the  moment  equations  of  motion  in  terms  of  the  Euler  angles  for  an  axisym- 
metric  vehicle  may  be  written^ 


M  «  10  +  I  piii  sin  6  -  Iip^  sin  0  cos  0  (2) 

n  x 

M  ■  I  ($  sin  G)  +  IGip  cos  0  -  I^pO 
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III.  AERODYNAMIC  MOMENTS 


It  is  assutaed  that  the  only  forces  acting  on  the  vehicle  that  contribute 


to  the  aerodynamic  moments  M,.,  M^,  are  normal  forces  from  angle  of  attack, 

#  4  «  1  4  ia  ^  mm  &  ak  a  *  J  hu4  m  ^4  t*  *  «u  4  u4  a  au  4  «•  »•  «  4  A  PA  4  •«  n  4  «««■  •*•«»> 

a  ▲at  aiU  awiuee  aavtu  Ctuii.au  opau  a  xuo  |  «Uu  pikuti  emu  /•*>  Uttwpau^  a5a  v.v»  • 

For  simplicity  only  two  spin  fins  are  considered,  although  the  results  so 
obtained  could  be  generalised  to  any  number.  The  fins  are  assumed  to  be  located 
on  the  aft  end  of  the  vehicle  a  distance  x  behind  the  vehicle  center  of  mass, 
as  shown  in  Fig.  2.  The  fins  are  canted  to  produce  a  positive  (clockwise) 
roll  moment  and  are  designated  Fin  1  and  Fin  2  in  order  to  differentiate  b.^tween 
the  windward  and  leeward  positions  at  any  instant.  The  roll  angle  $  describes 
the  roll  orientation  of  the  vehicle  relative  to  the  plane  of  total  angle  of 
attack  and  is  referenced  such  that  4>  -  0  when  Fin  1  is  leeward  and  Fin  2  is 
windward.  The  aerodynamic  moments  due  to  the  1  in  lift  forces  and  are  then 


M  -  <L.  +  L  )  ~ 

^f ins  1  1  1 


1  “  (L.  -  L„)  x  sin 

nfins  1 


■l  *  (L?  -  L.)  x  cos 

Cfins  1  1 


and  the  total  aerodynamic  moments  acting  on  the  vehicle  are 


Ms ' 


fi  -  M_ 


■^/  2C"  'I 

M  -  M  -  C  6qSx  -  l-C  +  — ~-|e 

n  ins  N0  6t  2u  \\  md2  / 


The  last  tens  In  the  second  and  third  of  Ena.  (4)  arm  nlceh  and  yaw  damn  1  no 
moments,  where  6  and  $  sin  6  are  the  pitch  and  yaw  rates ,  respectively. 

The  fin  lift  forces  are  given  by 


4  “  ca6(1)  Vs 

h  '  Ct.<2)  *2'* 

0 


(5) 


where  ^  and  ^  ^  are  fin  lift  derivatives  and  6^  and  6^  ere  the  effective 
6  £ 

cant  angles  of  FinB  1  and  2,  respectively.  The  effective  cant  angles  consist 


of  the  fixed  cant  angle  of  the  fins  relative  to  the  axis  of  symmetry  of 
the  vehicle  plus  the  component  of  angle  of  attack  that  may  add  to  or  subtract 


from  fiQ,  depending  on  the  roll  orientation  The  effective  cant  angle  can  be 
written 


6,  ■  6  -  6  sin  $ 

1  o 

5.  ■  6  +0  sin  $ 

/  o 


(6) 


The  dependence  of  the  fin  lift  forces  on  angle  of  attack,  is  expressed  through 
the  fin  lift  derivatives  according  to 


c,  (1)  -  C.  [1  +  Of.  (d>>] 

Ci  (2>  "  C2  11  +  0f2(<t>>1 
£  £ 


<7> 


¥ 


U  «»«)..>  --  J  •?  /  i\ 


where  C,  In  th«  94*Yft  Anal  n  nf 

- — « —  —  — - — -  *nu  l2vv^  *re  perxoaic 

functions  that  describe  the  windward  and  leeward  positions.  The  simplest 

first  order  dependence  on  angle  of  attack  is  assumed.  If  we  choose  for  ^($0 

and  £he  Enactions 


*l(4>)  ■  j  (a  -  cos  ♦) 


f?(4>)  •  -r  (a  +  cos  $) 


then  the  windward  and  leeward  values  of  C.  become 

<5 


<ci  >  -  c«  U  +  (^-—ue) 

6  windward  *6  ‘ 


<ci  >  "  C,  [1  -  (l-l-»)A0j 

6  leeward  *6  * 


(8) 


(9) 


For  a  ■  1,  the  windward  effectiveness  increases  with  angle  of  attack  while 
the  leeward  effectiveness  remains  unchanged,  whereas,  for  a  <  1,  the  windward 
effectiveness  increases  and  the  leeward  effectiveness  decreases  with  angle 
of  attack.  Equations  (5)-(8)  give,  for  the  fin  lift  moments  of  Eq.  (3),  the 
expressions 


M£  "  (1  +  Sin  <t>  cos  ^]C.  qSd 

fins  *5 
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(10) 


Mn  •  •In  ♦  eos  4 

fin* 

+  <2  *  aA8\  «<n2<ilr  «R*A 
*6 

H _  •  l<  A  cos2* 

Cfins  ° 

+  (2  +  «a6)  gin  *  cos  +)C.  qSxO 

l6 

For  angles  of  attack  greater  than  the  trla  values  from  configuration  and  aero¬ 
dynamic  symmetries  (assumed  to  be  sero  here) ,  the  windward-meridian  rotation 
rate  $  will,  in  general,  be  nonsero;^  l.e.,  the  vehicle  rolls  about  its  axis 
relative  to  the  wind.  Consequently,  the  sinusoidal  functions  of  $  in  Eq.  (10) 
can  be  averaged  over  a  cycle  to  yield  average  values  for  the  aerodynamic  moments 
due  to  the  unbalanced  fin  lift  forces*  Making  use  of  the  definite  Integrals 


1  r2*  1  r2ir  2 

•=—  F  sin  *  cos  <f>  d*  -  0,  ■=-  I  sin  *  d* 

211  Jq  2ir  •'o 


1  f2v  2  1 

*  2*  / Q  cos  ♦  **  2 


we  obtain  for  the  average  values  of  these  moments 


M, 


(1 


’fins 


6QqSd 


M  -  -(1  +  qSxe 

ins  2  l& 


(11) 


M. 


'fins 


“5  C,  6  qSx 
2  l6  0 


*» 
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IV.  EQUATIONS  OF  MOTION 


The  nomen t  aquations  of  notion  in  tsrns  of  the  Euler  angle  coordinates, 
obtained  from  subetitution  of  the  aerodynamic  moments  from  Eqs.  (4)  and  (11) 
into  Eq.  (2).  can  be  written 


(I  +  if*  •  i 


•%  aA  %  §  ••  •  *2 

•(I  +  — ~)t0  -  u»  0  -  v0  «  0  +.  ppip  sin  0  -  ^  sin  0  cos  0 


(12) 


4  6  TA0  -  vip  sin  0  ■  4r  Op  sin  0)  +  §<p  cos  0  -  yp0 
2  o  at 


where  the  coefficients  are 


ct64.’M 


,  c» 'Sx.t 

.2  a 


u> 


c#  qSx 
X6 


v  mT 


(13) 


2cn  1 

v-#:  [-c  +  — f- 

2Iu  1  "q  md2 


The  parameter  w  is  the  natural  pitch  frequency  of  the  vehicle. 
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V.  EQUIVALENT  MAGNUS  COEFFICIENT 
AND  YAW  DAMPING 

T^.6  o*h.»  M  Pn  n  1  \  Afiia ¥&A  Ko  Plw»  f <«  14ff  imKft  1  anr* 

•  ^  «k»  MtWMWM*  *  *  y  «M  ***f  *  <  **  r  y  W«  VW  •«««•  **  J  V  •  4i»  •••  «•  «MIW  *a*M»*V  V  f 

ins 

is  equivalent  to  a  Magnus  moment 

M  ■  C  p0  ( 

magnus  n  2u 

pa 


in  which  the  Magnus  coefficient  C  has  the  value 

npa 


C  -  AC.  6  (jf>  (~r) 
n  E6  o  d  pd 


Unlike  the  classical  Magnus  moment,  the  yaw  moment  M  is  independent  of 

**fins 

roll  rate  in  the  first-order  approximation.  The  fin  lift  derivative  C. 

6 

is  defined  such  that  the  product  C.  6  is  the  zero-angle-of -attack  value 

*6  0 

of  the  total  roll  moment  coefficient  induced  by  the  canted  spin  fins. 

0 

The  fin-induced  yaw  moment  can  also  be  related  r:o  an  equivalent  yaw  damping 
coefficient.  From  the  yaw  equation  of  motion  in  Eqs.  (12),  with  the  small-angle 
approximation  sin  0=0,  the  fin-lift  unbalance  term  is  equivalent  to  a  yaw 


damping  coefficient  °f  magnitude 


PRECEDIN6  PABE  BLANK 
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Following  the  approach  of  Kef.  6,  we  find  the  quasi-steady  value  of  the  precession 
rate  from  the  pitch  equation  of  notion.  Again  using  tint  snail  angle  approxi¬ 
mations  sin  C  =  S  and  cos  0  s  1,  we  can  write  this  equation  in  the  fora 

6  +  Iw2  +  (1  +  S—)  T  +  up*  -  *2)  6  v$  -  0  (17) 

which,  if  we  neglect  v  and  assuae  aA6/2  «  1,  yields,  for  the  two  precession 
nodes,  the  quasi-steady  values 

*♦.-  ■  “ !  V(*i)2  *  “2  * 1  <l8> 

The  term  T  under  the  radical,  which  arises  from  the  fin-lift  unbalance,  increases 
slightly  the  static  stability  of  the  vehicle.  Since  Eq.  (Id)  yields  one  positive 
and  one  negative  value  for  the  two  precession  modes,  the  equivalent  yaw  damping 
from  Eq.  (16)  will  be  either  positive  or  negative,  depending  on  which  precession 
mode  prevails.  The  positive  precession  mode  is  undamped,  whereas  the  negative 
precession  mode  is  damped,  which  is  consistent  with  the  results  of  Ref.  2. 
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VI.  REENTRY  VEHICLE  ANGLE-OF-ATTACK 
CONVERGENCE  ENVELOPE 

The  Influence  of  Che  fin  lift  unbalance  on  the  angle  of  attack  convergence 
of  a  rolling  reentry  vehicle  can  be  determined  directly  from  the  result  of 
Ref.  6,  with  the  equivalent  yaw  damping  coefficient  defined  by  Eqs. 

(16)  and  (18).  If  this  coefficient  is  added  to  the  usual  yew  damping  coefficient 
and  t  is  assumed  small  relative  to  ui  ,  we  obtain,  for  the  average  value  of 
the  angle-of-attack  convergence  envelope ,  the  expression 

0  2  I  1  f 

rm  d  +  0  >  exp  r  i  J0 

6  TA  1  ) 

-  j - -i  **  • 

Up  V  1  +  oz  J 

where  the  top  signs  in  the  t  and  f  signs  In  the  exponent  corresponds  to  the 

positive  precession  mode  and  Che  bottom  signs  correspond  to  the  negative 

precession  mode.  The  parameter  o  is  the  ratio  of  the  natural  pitch  frequency 

to  the  reduced  roll  rate  up/2.  For  a  slender  vehicle  with  u  »  1/10.  for  example 
2  /  2 

j  »  400/(p/w)  ,  which  is  a  strong  function  of  the  ratio  p/ui.  For  roll  rates 

2 

of  the  order  of  the  pitch  natural  frequency  or  less,  0  >>  1  and  Eq.  (19)  reduces 

to  the  simpler  form 

6_ 
u 

o 
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To  illustrate  the  effect  of  the  unbalanced  fin  lift  forces,  we  compare 
the  result*  ka,  (19)*  with  digital  computer  solutions  of  the  equations  of 
motion,  Eq.  (12).  Figure  3  shows  a  comparison  of  angle-of-attack  histories 
for  the  two  extreme  cases  of  exoatnospheric  motion  that  result  in  nearly 
circular  motion  in  both  the  positive  and  negative  precession  modes. ^  Also 
shown,  for  comparison,  is  the  engle-of-atcack  behavior  for  a  negative  pre¬ 
cession  case,  in  which  the  fin  lift  unbalance  is  nonexistent.  The  theoretical 
approximation,  Eq.  (19),  is  Included  for  comparison  in  each  case.  For  reentries 
in  which  the  axis  of  the  exoatnospheric  precession  cone  does  not  coincide  with 
the  velocity  vector,  the  initial  angle  of  attack  will  be  oscillatory.  The 
initial  average  value  of  the  oscillation  envelope  0^  for  such  cases  is  the 
angle  between  the  velocity  vector  and  the  axis  of  the  precession  cone. 

Figures  4  and  5  show  computer  simulations  of  two  such  cases  in  which  the 
velocity  vectors  lie  inside  and  tangent  to  the  exoatnospheric  precession 
cones,  respectively.  The  theoretical  approximations  to  the  average  values  of 
the  oscillation  envelopes,  computed  from  Eq.  (19),  are  shown  for  comparison. 
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Figure  3.  Comparison  of  Angle-of-Attack 
Histories 
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VII.  SIMMARY  AND  CONCLUSIONS 


A  Magnus- type  testability  or  tinned  Missiles,  caused  by  an  unbalance 
in  fin  lift  forces  from  unequal  effectiveness  of  windward  and  leeward  fins,  has 
been  described  quantitatively .  Although  of  different  origin  than  the  classical 
Magnus  weent  of  a  spinning  cylinder  in  a  cross  flow,  the  fin- lift-induced  yaw 
soaent  acts  in  the  seas  direction  and  causes  damping  of  negative  precession 
and  undamping  of  positive  precession  Motions.  Unlike  the  classical  Magnus 
aowsnt ,  the  fin-induced  yaw  Moment  is  Independent  of  roll  rate  in  the  first- 
order  linear  approximation.  The  equivalent  Magnus  coefficient  is  found  to 
have  the  value 


where  A  is  the  slope  of  the  fin  effectiveness  vs  angle-of-attack  curve,  assumed 

to  be  linear,  and  C.  is  the  fin-induced  roll  moment  coefficient  of  Che  missile 
*o 

A  closed-form  solution  has  been  obtained  for  the  angle-of-attack  convergence 
envelope  of  a  finned  reentry  vehicle.  This  result  indicates  the  magnitude  of 
the  fin  lift  unbalance  required  to  cause  s  net  undamplng  of  the  angle  of  attack 
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A  dynamic  instability  is  described  for  rolling,  finned  missiles  with 
canted  flna.  The  Instability  occurs  as  an  u?t damping  of  the  angle  of  attack 
and  is  caused  by  differential  lift  from  unequal  effectiveness  of  the  windward 
and  leeward  fins.  The  instability  is  similar  to  a  Magnus  instability  in  that 
a  yawing  moment  occurs  that  causes  damping  of  negative  precession  and 
undamping  of  positive  precession  motions.  The  equivalence  of  the  Magnus 
coefficient  and  the  differential  fin  lift  forces  is  derived,  and  the  insta¬ 
bility  is  demonstrated  with  computer  simulations  of  the  equations  of  motion 
for  the  angle-of-attack  convergence  of  a  reentry  vehicle.  The  computer 
results  are  compared  with  a  closed-form  solution  for  the  angle-of-attack 
convergence  envelope. 
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